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We study the motion of the Gaussian obstacle potential created by blue detuned laser beam 
through a phase-separated binary condensate in pancake-shaped traps. For the velocity of the 
obstacle above a critical velocity, we observe the generation of vortex dipoles in the outer component 
which can penetrate the inner component. This is equivalent to finite, although small, transport of 
outer component across the inner component. In the inner component, the same method can lead 
to the formation of coreless vortex dipoles. 



I. INTRODUCTION 

Generation of vortex dipoles in dilute Bose-Einstein 
condensates (BECs) when a superfluid moves past an 
obstacle with a velocity greater than a critical velocity 
(speed of sound locally) has been established by various 
numerical studies [I]-[S]. It has been established, both 
theoretically and numerically, that even the flow of the 
BEG below bulk sound speed can result in the flow ve- 
locity greater than the local sound speed at the equa- 
tor of the obstacle. It implies that the critical velocity 
for the creation of vortex-antivortex pairs is lower than 
the bulk sound speed. Above this critical speed vortex- 
antivortex pairs are created periodically, resulting in an 
oscillatory drag force. This is due to the phase differ- 
ence between the main flow and stationary wake behind 
the obstacle. Whenever this phase difference grows to 
2iT a vortex-antivortex pair is created. Recently, Neely 
et al. 1^ observed the generation of vortex dipoles in a 
pancake-shaped ®^Rb condensate, when the condensate 
was moved past a Gaussian obstacle, created by blue- 
detuned laser beam, above critical velocity. Real time 
dynamics of single vortex line and vortex dipoles has also 
been studied experimentally [7j . The flow of the miscible 
binary and spinor condensates across a Gaussian obsta- 
cle potential has also been investigated theoretically in 
Refs.[8 lOJ. Vortex dipoles or vortex rings are also pro- 
duced as one of the the decay products of dark solitons 
(a notch or density depression in the condensate with a 
phase slip of tt across it), due to the onset of a transverse 
instability in 2D and 3D geometries, namely snake insta- 
bility. The decay of dark soliton into vortex rings has 
already been observed experimentally [TTJ [12]. In this 
context, it should be noted that motion of the obstacle 
across a ID BEG generates grey solitons P^15j . The 
snake instability is also responsible for the decay of ring 
soliton into a necklace of vortex dipoles [TB] . 

In the present work, we study the flow of phase- 
separated binary condensates across a Gaussian obstacle 
numerically. We consider two experimentally realizable 
binary condensates: the first of these has ®^Rb and *^Rb 
and the second has the two hyperfine states of *^Rb as the 
two constituent species. In both these condensates, one 



of the scattering length can be tuned by using magnetic 
Feshbach resonances [T71 [TS] . This allows us to choose 
the tunable scattering lengths values suitable for gener- 
ation of coreless vortex dipoles. We observe the genera- 
tion of vortex dipoles for the flow velocity greater than 
the critical velocity. In our simulations while moving the 
obstacle, its strength is continuously decreased. We con- 
sider two cases: (a) the obstacle laser potential, initially 
in the outer component, is moved towards the center, and 
(b) the obstacle laser potential, initially in the inner com- 
ponent, moves towards the interface of the binary con- 
densate. With the successful realization of vortex dipoles 
in single component condensates using obstacle potential 
[B] , the present studies give more than enough reasons to 
extend these experimental studies to binary condensates. 



II. BINARY CONDENSATES AND VORTEX 
DIPOLES 

A. Meanfield description of binary condensates 

The dynamics of binary condensate at T = OK can be 
very well described by a set of coupled GP equations 



2m 



2 

V2 + X/,(r,t)+^t/.,|*,(r,t)p 



*»(r,t) 



(1) 

in mean field approximation, where i = 1, 2 is the species 
index. Here Ua = Airh^aa/mi, where m^ is the mass 
and an is the s-wave scattering length, is the intra- 
species interaction, Uij = 2'iTh?'aij /rriij, where m^- = 
mirrij / (rrii +mj) is the reduced mass and aij is the inter- 
species scattering length, is the inter-species interaction, 
and Vi{v) is the trapping potential experienced by zth 
species. In the present work, we consider binary con- 
densate consisting either of ^^Rb and ^^Rb for which 
mi ~ m2. Furthermore, we also consider identical trap- 
ping potential for both the species, the total potential is 
then ith species 
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where Vohs{.x,y,t) = Voit) exp{-2([x - xo(t)]^ + 2/^)/wo} 
is the potential created by blue detuned laser beam, and 
a and /3 are anisotropy parameters. To rewrite the cou- 
pled GP equations in scaled units, define the oscillator 
length of the trapping potential 



(3) 



and consider hu as the unit of energy. We then divide 
the Eq.Q by hu and apply the transformations 
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The transformed coupled GP equation in scaled units is 
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and 



where Uu 

= a^g^V^. From here on, we will represent the scaled 
quantities without tilde unless stated otherwise. 

In pancake-shaped (highly oblate) traps, trapping fre- 
quency along axial direction is much greater than the 
radial directions, i.e., /3 ^ a ~ 1; in which case Eq.([5]) 
can be reduced to two dimensional form by substitut- 
ing (/)(r,i) = ^p{x,y,t)^{z)exp{-i(3t/2) [19|. Here ^ = 
{(3/{27r)y^'^ exp(—/3 2:^/4), the ground state wave function 
in axial direction. The reduced two dimensional form of 
coupled GP equation in pan-cake shaped traps is 
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and Uij = 
work, we con- 
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sider ui2 > y^uiTuii^ so that the ground state of the 
binary condensate is phase-separated. The ground states 
of the binary condensates in the phase separated domain 
have been studied in Refs.pU]-}??] In pancake-shaped 
traps the binary condensates has cylindrical interface 
separating the two components |26] . 



B. Vortex dipole trajectory in a single component 
condensate 

The total velocity field experienced by each vortex in 
vortex dipole is the vector sum of two component fields. 
One of these fields is the velocity field generated on each 
vortex due to inhomogeneous nature of the condensate 
in the trapping potential. This inhomogeneity generated 
velocity field is responsible for the rotation of an off center 
vortex around the trap center [351 [5^ . In addition to this. 



each vortex creates a velocity field varying inversely with 
the distance from its center, which is experienced by the 
other vortex of vortex dipole. If (a;, y) and (x, —y) are the 
locations of the positively and negatively charged vortices 
of the vortex dipole, respectively, then the total velocity 
field experienced by the positively charged vortex of the 
vortex dipole 
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where Wpr is the rotational frequency of a vortex with 
charge -1-1 in the condensate. In terms of component 
velocities, the previous equation can be written as 
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This coupled set of differential equations describes the 
trajectory of the positively charged vortex of the vor- 
tex dipole. In the present work, we move the obstacle 
along X-axis, leading to the generation of vortex dipole 
located symmetrically about x-axis. Since y is small (and 
positive) for the positively charged vortex of the vortex 
dipole at the instant of generation, one can neglect —tUpj^y 
in first equation of coupled Eqsj8]to obtain 
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Eqs|9] imply that vortex initially to the left of the origin 
moves towards origin with increasing speed along positive 
X-axis but decreasing speed along negative y-axis. At x = 
0, dy/dt = and beyond this vortex moves away from 
origin with decreasing speed along x-axis but increasing 
speed along positive y-axis. Near the origin WprX <^ l/2y, 
then the velocity of the vortex can be approximated as 
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In this case the velocity of the vortex is mainly induced 
one, and it moves with induced velocity towards positive 
X direction. On the other hand if y is too large, one can 
neglect l/2y term in Eqjsjto obtain 
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In this case, the velocity of the vortex is mainly due to 
the inhomogeneous nature of the condensate, resulting 
in the rotation with frequency Wpr. The transition be- 
tween vortex induced velocity to inhomogeneity induced 
velocity occurs when dx/dt = 0, i.e., velocity along x-axis 
changes sign. From Eqsjsj it occurs at y = ^yl/2u!pl■. 
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III. OBSTACLE MODIFIED CONDENSATE 
DENSITY 

In the phase separated domain of binary condensates, 
the species with lower repulsion energy forms a core and 
the other species forms a shell around it. For convenience, 
we identify the former and later as the first and second 
species, respectively. With this labelling, interaction en- 
ergies uii < U22 and for equal populations this implies 
flu < a22- There is large difference of the local coherence 
length S^i — Xj \J2niUii when u\\ <^ U22- This results in 
different density modifications for the two species when 
an obstacle beam is applied. The dynamics is most sensi- 
tive to the coherence length difference while the obstacle 
is crossing the interface layer of the binary condensate. 

To analyze the density perturbations in presence of 
the obstacle beam, let i?in be the radius of the inner 
species or the interface boundary. And, let i?out be the 
radial extent of the outer species. In the absence of the 
obstacle beam, the chemical potential of first species is 
/ii = -Rfn/4 -I- uii/ [ttR^^). Here, the second term is cor- 
rection arising from the presence of the second species. 
The chemical potential of the second species is 
Both of these are in scaled units. We consider the ob- 
stacle beam initially {t = 0) located at (— i?out,0) and 
traversing towards the center with velocity fobs- As it 
moves, its intensity is ramped down at the rate rj. The 
location of the beam at a later time is 
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and intensity of the beam is 
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where Vo(0) is the initial intensity of the obstacle beam. 
At the starting point, the total potential V^(i?out, 0, 0) > 
i?Q^^/2 and the density of the outer species |V'2p is zero 
around the center of the obstacle beam as is shown in Fig. 
[ija). However, as it traverses the condensates with de- 
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FIG. 1. The stationary state wave functions of ^^'^Rb-^^Rb 
binary condensate obtained from imaginary time propagation. 
For (a) the obstacle is located near i?out, for (b) the obstacle 
is located at the interface, and for (c) it is well inside the bulk 
of the inner species. 

creasing intensity, at some later time t' , V{xo{t'), 0, 0) < 



■Rout/2 



Density |V'2p is then finite within the obstacle. 
For compact notations, hereafter we drop the explicit no- 
tation of time dependence while writing Xo{t) and Vo{t). 



A. Density modification at interface 

As the obstacle beam approaches the interface, the in- 
tensity is weak to expel out the outer species. However, 
when the beam enters the inner species with any finite 
intensity V{Rin,0,t) > fii, it is sufficient to make IV'iP 
zero around the beam [see Fig. [ija)]. The interface is 
modified and equipotential curve of the interface is 
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A round the center of the obstacle beam, 
\/ix — xq)'^ + < w/V2, and the obstacle poten- 
tial is 
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Hence from Eq. ( 14 ) the equipotential curve in the neigh- 
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Equipotential curve in the vicinity of the obstacle is thus 
modified into a circle with center at (— C/2A, 0) and ra- 
dius i?in. In our numerical calculations, the Gaussian 
obstacle potential is quite narrow and u'2/(4Vb) ^ 1. 
Retaining only the leading order terms in /AVq, the 
center is approximately located at [xq -{- uP'Xf)/ {AVq),Q) 
and radius is 




(17) 



Sufficiently away from the obstacle \J {x — xqY + J/^ > 
110/^/2 and the obstacle potential is negligible. The 
equipotential curve then remains unchanged, — x"^ + 
. Thus an approximate piecewise expression of the 
equipotential curve is 
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This curve defines the interface geometry of the inner 
species in the TF approximation. Surface effects, how- 
ever, smoothen the density profile around the meeting 
point of the two regions as is shown in Fig. [l]^b). The 
obstacle beam is repulsive and expels the condensate 
around it such that bulk density is regained over a dis- 
tance of one healing length. 



B. Fully immersed obstacle beam 

A critical requirement to form coreless vortices is com- 
plete immersion of the obstacle beam within ni. Based 
on the previous discussions, as the beam approaches the 
origin, the last point of contact between the beam and 
interface at Ri^ lies along x-axis. To determine the condi- 
tion when complete immersion occurs, consider the total 
potential along x-axis 
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where, the Gaussian beam potential is considered upto 
the second order term. The expression is appropriate in 
the neighborhood of the beam and has one local minima 
and maxima each. There is also a global minima, how- 
ever, it is not the correct solution as it lies in the domain 
where ^ > Correct global minima is the one located 
at X = and associated with the harmonic potential. 
The obstacle is considered well immersed when the local 
minima is located at the interfacial radius Rin- The local 
minima is a root of the polynomial 



x^ + a2X^ + aix + oo = 0, 
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obtained from the condition dV/dx = 0. Here, 
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For the trapping potential parameters considered in the 
present work, all the three roots are real. The root cor- 
responding to local minima is 
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and Xmin is a function of time as the obstacle potential 
is time dependent. However, for compact notations, like 
in Vq and xq, hereafter we drop the explicit notation of 
time dependence while writing x^in(t). Value of Xmin is 
real when it satisfies certain conditions |30| and these are 
met for the experimentally realizable parameters. The 
obstacle is then completely immersed when x^i^ — 
and let denote the time at which this condition is 
satisfied. Once the obstacle beam is well inside the inner 
species, within the beam ni is zero but n2 is nonzero. 
It then forms a second interface layer, assisted by the 
obstacle beam, which embeds a bubble of the 712 within 
Til. Recollect, the first interface layer is located at Rin 
and it is where ^2 encloses ni. The second interface, 
unlike the one at is a deformed-ellipse and label it 
as F. It touches x-axis at two points and Xmin is the one 
farther from the origin. From TF-approximation, the 
density of n2 within F is 
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where the subscript F is to indicate the density is within 
F. At ti„i, n2r is in equilibrium with the bulk of 712, and 
let the number of atoms in the be N2r, which can be 
calculated from n2r- 

Except for deviations arising from mean field energy, 
the second interface F closely follows an equipotential 
surface. From the conditions of pressure balance at two 
interfaces 
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The first equation is valid all along the interface at R^n- 
However, the second is derived at (rniin,0), but to a 
good approximation the condition is valid along F as it 
is close to an equipotential surface. The density of the 
first species at (xminjO) is 
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Similarly density of second species at at (x,„in,0) is 
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Inside F, where the obstacle potential is dominant, rii 
decays exponentially. Along the x-axis 
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where, (5x_ is the distance of the point measured from F. 
Similarly the density of the second species just outside F 
is 



n2((5x+,0) = n2(xniin,0)exp ( --^6x4 
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where 5x+ is distance from T but away from the obstacle 
beam, and A2 is the penetration depth of the second 
species. In these expressions the penetration depth is 
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Based on the densities at the interface region, we can 
incorporate the interspecies interaction and calculate the 
effective potential of each species. 



C. Obstacle assisted bubble 

At a time At after the obstacle is immersed in rii, the 
location and amplitude of the obstacle potential are 

2;o(tim + At) = -Rout + l'ob(iim + At), (28) 
yo(tim + At) = Vo{0) + Vobitin, + At). (29) 

Equilibrium TF n2 within the obstacle potential at this 
instant of time is 
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This, however, is higher than the density distribution at 
tim, that is n2r{x,y,tini + At) > n2r{x,y,t^^) as the 
potential V is lower. This is on account of two factors: 
first, the amplitude of the obstacle potential decreases 
with time; and second, the harmonic oscillator potential 
is lower at a;o(^im + ^t)- The number of atoms, however, 
does not change from the value at ti^ unless there is a 
strong Josephson current. Density n2r is thus below the 
equilibrium value once the obstacle beam is well within 
rii . This creates a stable bubble of 712 assisted or trapped 
within the beam and is transported through the rii. 

Departure of n2 from the equilibrium is not the only 
density evolution within the beam. There is a progressive 
change of nir as the beam moves deeper into ni. At 
time tini, when the obstacle is completely immersed in ni 
the effective potential, experienced by ni, V{x,y,tiin) + 
n2TUi2 is larger than /ii. So, ni is zero within the beam. 
However, if the rate of ramping 77 is such that at a later 
time V(x,y,ti^ + At) + n2rui2 < Mii while the beam is 
still within rti, there is a finite ni within the beam. 

/ii - y, + At) - Mi2n2r ,„,s 
nir = (31) 

Ull 

As the growth of nir is in its nascent stage, Eq. (pO|)is 
still applicable. Since ai2 > -^011022 for the condensate, 
in TF approximation both uiy and n2r can not be simul- 
taneously non-zero. At the same time, 7i2r is forbidden 
to migrate to the bulk n2 due to the ni generated poten- 
tial barrier in the region between interfaces F and i?in. 
To accommodate both rii and n2 within the beam, the 
shape of interface F is transformed to increase n2r- So 
that 712 is zero in certain regions within the beam where 
the condition V{x,'y,tnu + At) + 7i2rui2 < Mi is satis- 
fied. This mechanism is responsible for obstacle assisted 
transport of 77.2 across ni. 



D. Obstacle induced density instability 

The harmonic potential increases once the obstacle 
crosses origin, but the obstacle potential continues to de- 
crease. A fine balance is achieved if 
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the dynamical rate of change of obstacle potential, is 
same in magnitude and opposite to the rate of change 
of the harmonic potential dV^o/dx. It is, however, dif- 
ficult to satisfy this condition across the entire obstacle 
potential as the Gaussian beam has has different profile 
than the oscillator trapping potential. But, the param- 
eters may be tuned so that the condition is satisfied at 
the center of the beam, then 
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In the present case 7; is a constant and the condition in 
Eq. ( 33 1 is satisfied only when the beam is at a specific 
point Xe, then 
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Before reaching Xc, the potential within the beam de- 
creases, but increases after crossing it. At some point 
beyond Xc, the total potential experienced by 712 bub- 
ble can become greater than fi2, resulting in an unstable 
bubble. 



IV. STABILITY OF CORELESS VORTEX 
DIPOLE 

To analyze the energetic stability of coreless vortex 
dipoles in phase-separated binary condensates, we adopt 
the following ansatz for the wave function of inner and 
outer species 



ipi{x,y) = a[i(x- vi) + (y^ - vl)] exp ^ 
■02 (a;, y) = 6 (x^ -f y^) exp f ^ 
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where {vi, ±^2) is the location of the vortex dipole. With 
this choice of ansatz, the vortex dipole in the inner species 
is a normal vortex dipole. For simplicity of analysis, con- 
sider the trapping potential is purely harmonic, i.e. 
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Normalization yields the following constraints on a and 
b 
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and reduces the number of unknowns to two vi and W2- 
The energy of the binary condensate using the ansatz in 
Eq. (lasl) is 
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+4vf - + 4v|) + \?ux2 (15 + 8v\ - I2vl + 8w|)] 
+a2un [177 - 304^2 + 2^Qv\ + 416i;^ - 256w^ 
+256w^ + 32w^l5-8w^ + 16t;^)]}. (38) 

In the same way, we can also determine the chemical 
potentials of the two species in terms of the same param- 
eters. The stationary points of the energy function are 
the solutions of the coupled equations dE/dvi = and 
dE/dv2 = 0. These equations are solvable analytically, 
but the solutions are lengthy and cumbersome. The so- 
lution of particular interest is the vortex dipole located 
along the y-axis {vi = 0). This is consistent with the 
observation that, regardless of the orientation, from the 
symmetry the axis passing through vortex dipole can be 
considered as the y-axis. A rotation transformation is 
sufficient to achieve this. For the case considered, the 
energy extrema conditions simplify to equations 
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For the sake of illustration, consider, A^i — 10^, N2 — 
Ni/2, an = 0.99ao, 022 = 4.6ao, ai2 = 2.14ao, /3 ^ 
11.25, w/27r = 8, and Vo{t) = 0. Call this set of param- 
eters as parameters set b. The stationary points corre- 
sponding to minimum energy for this parameters set are 
{vi,v2) = (0, ±0.9678). 



A. Coreless vortex dipole 

To compare the energy of this normal vortex dipole 
with that coreless vortex dipole, consider that the cores 
of the vortices at (wi,W2) and (wi,— W2) are filled with 
outer species. We approximate the cores as circular with 
a radius of ^ = Xj ^2[i\^ the coherence length of the inner 



species. We approximate the density profiles of the outer 
species within the cores with the TF distribution 
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where Tc = ^/{x — vi)"^ + {y ± ^2)^ is the radial dis- 
tance from the core of each vortex and antivortex. The 
above expression for 'ip2 is renormalized to obtain the 
modified value of b, which in turn enables us to cal- 
culate the energy for the coreless vortex dipole. To 
obtain the analytic expression for energy, we integrate 

b (a;^ -I- y^) exp ^— ^ "^^ ^ over the whole real space with- 
out excluding the core regions. The modified value of b 
after the normalization is 
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With this approximation, the number of atoms of outer 
species embedded within the cores is 

bN2 = A^2 ^-^^^ — ■ (42) 

2^22 

These contribute to the total energy through the poten- 
tial and intra-species interaction energies. The energy 
contribution from the cores is then 

^ + -ivl + %evl + 3t'2' + 6«2 + ^;2) 
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The total energy of the binary condensate with coreless 
vortex is 



Er 
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where E is defined by Eqs. (38) and (41 1. For the pa- 
rameter set b, the energy per particle |i?c/(A'^i + ^^2)] 
is 4.18 which is smaller than the value of 4.26 obtained 
from Eqs. (37) and ((SSl). This implies that coreless vortex 



dipole is energetically more favorable than normal vortex 
dipole with the percentage energy difference of 1.9%. 

The Gaussian form of the ansatz imposes limitations 
on the domain of applicability. But, the ansatz pro- 
vide qualitative information on the stability coreless vor- 
tex dipoles. To relate to experimental realizations, con- 
sider ^^Rb-*^Rb binary condensate with an — 460ao, 
^22 = 99ao, and 012 — 214ao as the scattering length 
values and 2Ni ~ N2 — 10^ as the number of atoms. 
Here an is tunable with magnetic Feshbach resonance 
|31) . With this set of parameters, the stationary state 
of *^Rb-®''Rb binary condensate is just phase-separated. 
The trapping potential and obstacle laser potential pa- 
rameters are same as those considered in Ref.l6 , that 
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is w/(27r) = 8Hz, a = 1, /? = 11.25, Vo(0) = 93.0huj, 
and wq = lOfiui. Hereafter we term this set of scatter- 
ing lengths, number of atoms, and trapping potentials as 
set a. For this parameter set, in absence of the obsta- 
cle potential, the stationary point of minimum energy is 
(^'1,^2) — (0, ±0.75275). And, the difference in the en- 
ergies per particle of the binary condensate with normal 
and coreless vortex dipoles is 0.11. 



The energy can be minimized under the constraint that 
^1 and fj,2 satisfy Eqs. ( 47p8 ) to completely determine 
the stationary state of binary condensate with a normal 
vortex at the origin. 



2. Coreless Vortex 



B. TF approximation single vortex 



For a coreless vortex, the wave functions of the two 
species in TF approximation are 



As mentioned earlier, the applicability of the varia- 
tional antsatz is limited to the weakly interacting do- 
main. Experiments, however, mostly explore domains 
other than weakly interacting. To analyze stability in 
this domain, we resort to the TF approximation. 



1. Normal vortex 

To simplify the analysis, consider a normal vortex lo- 
cated at the center of the inner species. For a normal 
vortex, the wave functions of the first or inner species in 
TF approximation is 



l^ii-Vir 



Uii 



(45) 



Similarly, the wave function of the second or outer species 
is 



■02(7-) 



/m 2 - V{r) 
U22 



Rt <r'<R. 



> R, 



out 
2 

out 



(46) 



here V{r), i?out, and ^ are the trapping potential, 
radius of the interface between inner and outer species, 
radial extent of outer species, and coherence length of 
inner species, respectively. The chemical potentials ^1 
and IJL2 are parameters of the calculations and obtained 
from inverting the normalization condition 



N2 = 



x(1 + 4m?)] 



4u 



[2Rl (/.I -f - 2Rtlil + In (1 + 2Rlfii) 

(47) 
(48) 



22 



The total energy of the vortex with a normal vortex lo- 
cated at origin is 



E : 



■01 (r) 



r dr dr^ 2 



Mr) + Mr) [V{r) + ^|02(r)p] Mr)} , (49) 



ll^i-Vir 



■01 (r) 



^ Rl 



> Rl 



(50) 



Mr) 



I fi2 - V{r) - ui2'0i(r)2 



"22 



r'^e 



ll^2-V{r) 
U22 



Rt ^r'^R: 



> RL 



i,l) 



Normalizing the wave functions gives the constraint equa- 
tions in terms of fii and fj.2 as 

N, = [2i?f„ (Ail + 4^*?) - 2i?f„^2 + In (1 + 2R^^fi,) 

Swii^i 



No 



I6U11U22M1 

?4 ,,2 



{"21 [-l + ln4 + 8(ln2-l)/i2] + 



{4Ri^il + 8^1 M2 - l&RlAt^2 + 16Ai?M2 - 1) } (52) 



From the expressions, we can calculate the total energy 
of the system. It is a complicated expression which we 
refrain from writing. For the present case, it is perhaps 
sufficient to say, it is possible to obtain an analytical but 
cumbersome expression of total energy. 

For the parameter set a with Vo(t) = 0, binary con- 
densate with the coreless vortex has lower energy than 
the one with normal vortex. This is evident from Figj2] 
where the variation of total energy is shown as a function 
of i?i„. However, the stability vanishes if one considers 
Ni — N2 — 10^, ail = 51ao, 022 — 99ao in set a. It 
is then the normal vortex which is energetically stable. 
This is evident form Figj3| These results are in very good 
agreement with numerical results. This TF analysis, al- 
beit, is for single vortex clearly illustrates that depend- 
ing upon the interaction parameters, a phase-separated 
binary condensate can either support coreless or normal 
vortex. It must be emphasized, the energy, chemical po- 
tentials, and interface radius obtained from the TF ap- 
proximation are in very good agreement with numerical 
solution of coupled GP equations. For example, the nu- 
merical value of energy 66.32 corresponding to Fig|3] is 
very close to the value 66.06 from TF approximation. 
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85.52 



85.48 



85.46 



85.44 



Normal vortex 
Coreless vortex 




9.3 9.4 

Rin (units of aosc) 



9.5 



9.6 



FIG. 2. The energy of the binary condensate with Vo ~ and 
rest of the parameters same as those in the parameters set 
a as a function of Rin- The condensate has a vortex dipole 
located at origin. Black and blue curves are for coreless and 
normal vortex dipoles respectively. 




Rin (units of Eosc) 

FIG. 3. The energy of the binary condensate with A'^i = 
N2 = 10^, ail = 51ao, 022 = 99ao, Vo — 0, and rest of the 
parameters same as those in parameters set a as a function of 
Rin- The condensate has a vortex located at the origin. Black 
and blue curves are for coreless and normal vortex dipoles 
respectively. 



C. TF approximation with vortex dipole 

We use the TF approximation to examine the stabiUty 
of coreless vortices in the strongly interacting domain. 



1. Normal vortex dipole 

Assume the vortex affects the density of the condensate 
only within the core regions. We then adopt the following 
ansatz for binary condensate with a normal vortex dipole 



at (?;i,±t;2)- 



V'i(r) 



r 





= < 



■02 W 



//ii - V{x,y) 

Ull 

I /i2 - V{x, y) 
U22 



[{x~v,f + {y±V2r]^e^ ^ 
r , 9 9 53 

\ [ix-v,r+{y±v2r]>e 
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Rl^ix^+y^)^Rl, 

{x^ + y^) > Rl,, 
(x^ + y^) < Rl. 



(54) 



The vortex dipole contributes mainly through the kinetic 
energy of -01 , which may be approximated with the value 
of single species condensate given in Ref . [32] 



Vd 



2/i 



2V2 



(55) 



Using these ansatz the number of atoms are 

^ (1 + 4«2^i + Av^fi, - 8/i? - IRiA + 8i?2„Mi) 



N2 = 



8miiMi 



TT {Rl - 2/i2) ■ 



(56) 



In a similar way, we can evaluate the energy of the entire 
condensate. 



2. Coreless vortex dipole 



For coreless vortex dipole, we adopt the ansatz 



■01 (r) 



r 





x^+y' > Ri 

[ix-v,f + {y±V2f]^e 



lfii-V{x,y) \x^+y^^Rlk 



(57) 



Ull 



■02(f) 



I H2 - Vix,y) 
U22 



[ [{x-v,f + {y±V2f]>e, 
( Rl ^ (x' + y') ^ Rl^, II 

\[{x-v,r+{y±v2r]^e 

{x' + y')>RL, (58) 

{ x^ +y^ < Rl & 
\ [{x-v,f + {y±V2f]>e- 
Using these ansatz the modified expressions for N2 is 

N2 = TT^^ (2i?f„/i? + ^^il^Ji2 - ^Rlt4i^2 + 8/i?M2 - 1 

(59) 



8w22Aii 

-^vlni - 4w|^i) 



We can also calculate the total energy E of the system. 
The important change in E is the inclusion of interface 
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interaction energy Ei^t ■ It arises from the interface inter- 
actions at the cores of the vortex and antivortex. Based 
on Ref. [H], 

i^int = IPbTT^ - l) ' (60) 

3 V^aiia22 / 

where P is the pressure on the circumference of the cores 
and 

6^2 [ 3(mi +M2)7^ 1^/^^ 
4^1^2(012 - ^/alla22)J 

As in the previous section, energy can be minimized with 
the constraint of the fixed number of atoms. For the 
parameters set a without obstacle potential, the coreless 
vortex dipole has lower energy than the normal vortex 
dipole and is shown in Fig.|4]for the vortex dipole located 
at (0,±1). 

85.45 I ^ : 1 




9.2 9.4 9.6 9.8 

Rin (units of ao3c) 



FIG. 4. The energy of the binary condensate with Vo ~ and 
rest of the parameters same as those in the parameters set 
a as a function of 7?in. The condensate has a vortex dipole 
located at (0,±1). Black and blue curves are for coreless and 
normal vortex dipoles respectively. 



For Ni = N2 = 10^, flu — 51ao, a22 = 99ao and rest of 
the parameters same as in parameter set a, condensate 
with the normal vortex dipole has lower energy than the 
one with coreless vortex dipole (see Fig. [5]). These results 
are again in good agreement with the numerical results. 

V. NUMERICAL RESULTS 

For the numerical calculations, we again resort to ex- 
perimentally realizable system and paramteres. We again 
work with the ^^Rb-*^Rb mixture and parameter set a 
introduced earlier. Recall, the values of paramteres in 
this set are an = 460ao, 022 = 99ao, and 012 = 214ao 
2Ni = N2 = 10^ w/(27r) = 8Hz, a = 1, (3 = 11.25, 



66.6 




Rin (units of ao3c) 

FIG. 5. The energy of the binary condensate with A'^i = N2 = 
10**, ail ~ 51ao, 022 = 99ao, Vo = 0, and rest of the param- 
eters same as those in parameters set a as a function of Rin- 
The condensate has vortex dipole located at (0,±1). Black 
and blue curves are for coreless and normal vortex dipoles 
respectively. 

Vo(0) — 93.0/iw, and Wq — lOfim. As mentioned earlier, 
flu is tunable with magnetic Feshback resonance |31| . 

The transport of ^2 within the region of influence of 
Vohs is crucial to generate coreless vortices. As V^bs en- 
ters ni, it nudges ni and creates an interstice, which is 
occupied by 712- This is an outcome of the difference in 
coherence lengths (^^ = l/^/2niUii). The smaller ^ of 
n2 ensures that it recovers bulk value over shorter dis- 
tances. It implies that it is comparatively more difficult 
to displace the n2 than ni. As the obstacle reaches the 
interface with diminished strength, it is unable to dis- 
place n2, but not for ni due larger ^. This results in an 
inward protrusion of the interface, which detaches from 
the interface at a later time. For detachment to occur, 
it should be energetically favorable. The difference in 
f coupled with 012 ^ ^^0102 ensures that the interstitial 
filling of rt2 is energetically lower, in other words energet- 
ically favorable. This is easily confirmed with imaginary 
time propagation of GP equation to obtain the station- 
ary state geometry of the binary condensate with Vohs 
located inside ni as is shown in Fig. [ijc) . 

A. Motion of V^bs towards interface 

We first consider the motion of Vohs, initially located 
in the n2 towards the interface. For this case, we con- 
sider ^^Rb-®^Rb binary condensate with parameter set a 
and Vohs of Vb(0) = 125.0. The Vohs is initially located 
at a; = — 15aosc and moves with the speed of 180/im/s. 
It progressively decreases in strength and vanishes at 
X = Soosc- The obstacle creates a normal vortex dipole 
as it traverses 712- And, as Vobs approaches the interface, 
it transports the vortex dipole. Further motion of Vobs 
in the ni generates coreless vortices in the inner com- 
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FIG. 6. Coreless vortex dipole shedding when the obstacle 
potential traverses from the edge to the inner component in a 
binary condensate with Vb = 125, a;o(0) = — 15aosc and rest of 
the parameters same as those in set a. Figures a, b and c are 
the absolute values of wave functions at 0.21s, 0.33s, 0.56s, 
respectively; while figures d, e and f are the corresponding 
phase plots. 



ponent. The formation of the coreless vortex dipole is 
shown in Fig. [6j The plots in the figure show both the 
density and phase distribution. From the phase plot, it 
is clear that the phase singularity of the vortex is asso- 
ciated with Tlx- And, from the density plot it is evident 
that the vortex-antivortex cores are filled with ni. One 
effect which precedes the seeding of coreless vortex is the 
indentation of the inner interface after it cleaves when 
the 14bs is fully immersed in n\. This may be explained 
in terms of the dynamical evolution around T4bs and an- 
alyzed in two different ways. 



1. Inertia 0/712 

Just after the interface cleaves, when the Vobs is set 
into motion at t = tg, the n-i due to inertia, does not 
follow T4bs initially. As a result at to + neglecting 
the back fiow of 712, outer species (along x-axis) at x' G 
[2;miii(to), a;min(to + (5^)] is approximately same 



(62) 



here bt is time taken by the outer species to redistribute 
itself according to the displaced obstacle potential. As a 
result over a length a;min(to + <5i) — a;min(to), the density 
of inner species in TF approximation is 

ni{x ,0,to + ot) — . (63) 



Hence inner ni flows into the above domain and replaces 
n2. The circular symmetry of V^bs coupled with mo- 
tion along X-axis ensures that the aforementioned pro- 
cess starts along x-axis earlier; also the replacement of 
n-2 by ni takes place at the much faster speed as com- 
pared to the speed of the obstacle. This leads to creation 
of dimple on the inner interface. 



2. Pressure balance 

Let xo and T^^ be the location of the obstacle just after 
the interface cleaves and approximately circular inner in- 
terface centered around xq with radius Rm, respectively. 
After time St, let Xq = Xo+Sxq and F^.^ be the location of 
the obstacle and inner interface, respectively. Assuming 
that Tj.'^ has same radius as Txg, i.e. Rin, the points of 
intersection of Tm and F^' are 



5x1 



Thus the region ST defined as 

{x - Xof + < i?2^ x'of + 

with area 

Sxo 



(64) 



(65) 



A = nRt 



2R:„ COS 



2Ru 



Sxq 
2 



4i?f, 



Sxl, (66) 



lies within T^g but outside F^.' . Let Pi = n\uii/2 and 
P2 = 712^22/2 be the pressures exerted by inner and outer 
species on the left arc of 5T region, respectively. The po- 
tential experienced by ni along the left arc of (5F, which 
is approximately equal to harmonic trapping potential, 
does not change as the obstacle is moving away from this 
region, and so does the Pi. The potential experienced 
by n2 over (5F, which is equal to sum total of harmonic 
trapping and obstacle potentials, decreases as the obsta- 
cle moves from xq to xq -t- 5xq. In fact, as one moves 
from (a;',0), where x' e {xq — i?in,a;o -t- ^a;o — -Rin), to 
the (xcUc) along a curve with same sign of curvature 
as bounding arcs of ST region, the potential experienced 
by 712 decreases. This leads to the redistribution of 7i2, 
i.e. it moves to the regions with weaker trapping poten- 
tial. The redistribution will lead to 6P = Pi - P2 > 
which decreases as one moves along the aforementioned 
curve. This will lead to indentation of the inner inter- 
face. In our simulations, when the obstacle is moved, 
along X-axis, from a point near the center towards the 
outer interface, indentation is suppressed. It confirms 
that this mechanism indeed aids the indentation of the 
interface. 



B. Obstacle motion towards edge 

In this case, we start with the stationary state solution 
of ^^Rb-^^Rb binary condensate for the parameters set 



11 



a with obstacle laser potential maxima located at a; = 
— S.Oflosc- The static solution in this case is shown in 
first column of FigjTj which has roughly doughnut shaped 
density distribution of the outer species along with non- 
zero density in the circle of influence of the obstacle. 
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FIG. 9. The generation of normal vortex dipoles in the *^Rb- 
**^Rb binary condensate with an — 51ao, A''i = ~ 10^, and 
the rest of the parameters same as those defined in parameters 
set a. The obstacle potential, initially located at a:: = Oflosc, 
is moved with a velocity of 200/im/s up to a; = Soosc- First, 
second, and third columns are the absolute values of wave 
functions at t = Os, t = 0.01s, and t = 0.14s, respectively. 



FIG. 7. The generation of triply charged vortex dipole in 
the *^Rb-*^Rb binary condensate with parameters set a. The 
obstacle potential, initially located at x = — S.Ofloac, is moved 
with a velocity of 220/im/s up to s = Sflosc- First, second, 
and third columns are the solutions at t = Os, t = 0.14s and 
t — 0.32s, respectively. 



Such density distributions are prerequisite to create 
coreless vortices in the inner species. The obstacle 
is moved with a uniform speed of 220/zm/s, progres- 
sively decreases in strength, and finally becomes zero at 
X = S.Oaosc- With this velocity, we observe the creation 
of triply charged coreless vortex dipole as is shown in 
last column of FigjTj Singly and doubly charged vortex 
dipoles can be generated by moving the obstacle poten- 
tial at lower velocities. For example, FigjSj shows the 
singly charged coreless vortex dipole and the trajectory 
followed by it for the ^^Rb-*^Rb binary condensate with 
parameters set a. The obstacle laser potential moves 
from X — Oflosc Bi t = Os to X = S.Ooosc with the speed 
of 200Aim/s. 



We also considered the ^^Rb-® Rb with an = Slap, 
Ni — N2 — 10^, and the rest of the parameters same 
as those in parameter set a. The obstacle laser potential 
is initially located at origin x — Ooosc- The stationary 
solution in this case is strongly phase separated (ai2 S> 
^0102), and the circular region of influence of obstacle 
potential is devoid of atoms of outer component as is 
shown in first column of Fig|9j In our simulation, the 
obstacle potential potential traverses a distance of 5aosc 
along x-axis with a speed of 200/xm/s. We observe that 
the vortex dipole generated in the inner component is a 
normal vortex dipole with empty core. This is consistent 
with the TF results shown in Fig. j5j 



VI. CONCLUSIONS 




FIG. 8. The figure on the left side shows the vortex dipole 
in ^^Rb-^'^Rb condensate with parameters set a a,t t — 0.12s. 
The white curve indicates the trajectory of lower vortex. The 
figure on the right side shows the variation in position of the 
upper vortex with time. 



We have studied the motion of the Gaussian obstacle 
potential through a phase-separated binary condensate 
in pan-cake shaped traps. The motion of the obstacle 
leads to the generation of either coreless or normal vor- 
tex dipole depending upon the interaction parameters. 
The interaction parameters suitable for coreless vortex 
dipole can also be employed to have obstacle assisted 
transport of one species across another. We have used 
both a variational ansatz based method, applicable when 
Niaa/aosc ~ 1, and TF approximation based method, 
applicable when Nitta /ttosc 3> 1, to analyze the energetic 
stability of coreless versus normal vortices. Our stud- 
ies show that using Feshback resonances to tune one of 
the scattering length, e.g. scattering length of ^^Rb in 
^^Rb-^^Rb, it should be experimentally possible to create 
coreless vortices in phase-separated binary condensates. 
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